Dark matter effects in vacuum spacetime 
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We analyze a toy model describing an empty spacetime in which the motion of a test mass (and 
the trajectories of photons) evidence the presence of a continuous and homogeneous distribution of 
matter; however, since the energy-momentum tensor vanishes, no real matter or energy distribution 
is present at all. Thus, a hypothetical observer will conclude that he is immersed in some sort of dark 
matter, even though he has no chance to directly detect it. This suggests yet another possibility of 
explaining the elusive dark matter as a purely dynamical effect due to the curvature of spacetime. 
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I. INTRODUCTION 

According to the so-called standard cosniological 
model, the dustlike matter constitutes approximately 
27% of the total energy density of the Universe (the 
remaining 73% is dark energy, driving cosmic accelera- 
tion). All astrophysical and cosmological evidence in- 
dicates that more than 80% of the matter content of 
the Universe is in the form of some sort of nonluminous 
component (dark matter) [1, 2]. However, after almost 
80 years since the first indications appeared of the exis- 
tence of such a component [3, 4], the true nature of dark 
matter is still shrouded in deep mystery. Some fraction 
of the missing mass may be composed by baryonic ob- 
jects which do not shine, such as neutron stars, black 
holes, white dwarfs, brown dwarfs, and planets. How- 
ever, data from big bang nucleosynthesis and from the 
cosmic microwave background indicate that the overall 
baryon contribution, both luminous and nonluminous, 
cannot exceed about 20% of the total [5-7], which leaves 
the remaining 80% unexplained. 

The most credited hypothesis, suggested by the open 
questions raised by the standard model of particle physics 
and the resulting need of its extension, is that the non 
baryonic dark matter is composed by stable, neutral, cold 
heavy particles which interact only very weakly with or- 
dinary matter, such as several sorts of weakly interacting 
massive particles (WIMPs), axions, sterile neutrinos, and 
others. See [1] for an extensive review of the expected 
properties and methods of detection of such various pos- 
sible dark matter candidates. 

Several experiments and observations are under way in 
laboratories all over the world aimed at detecting signs of 
the existence of these elusive particles. See, e.g., [2] for an 
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updated review of the present state-of-the-art direct and 
indirect detection methods for dark matter. However, in 
spite of all these efforts, no conclusive evidence has so 
far turned up, though there are a few scanty signals that 
may be tentatively interpreted as dark matter signatures. 

In support of the hypothesis of the existence of stable 
WIMPs as being responsible for dark matter are theoret- 
ical motivations, such as the need for a solution to the 
gauge hierarchy problem [1], according to which these 
particles should have masses in the weak scale range 
(around 10 GeV - 1 TeV). Indeed, such a mass scale is 
just the right one for such particles to be produced with a 
relic density consistent with that required to explain the 
observed dark matter abundance and distribution (the 
"WIMP miracle"). 

There is also much hope that, if WIMPs are indeed 
there, they may soon be produced in pp high energy in- 
teractions at the LHC, thus turning up as missing mass 
signals in the shower of particles produced in a collision. 

An alternative approach to the dark matter problem 
is based on the modified Newtonian dynamics hypothe- 
sis (MOND). These models attempt to solve the miss- 
ing mass conundrum by modifying Newton's gravita- 
tional law at large distances, without introducing new 
forms of matter [8] . One of the relativistic realizations of 
the MOND paradigm is the tensor-vector-scalar gravity 
(TeVeS). This theory modifies gravity using a set of ad- 
ditional interacting fields to mimic dark matter [9] . How- 
ever, TeVeS has been recently ruled out by weak lensing 
observations, even though other MOND models are still 
permitted [10]. 

Another approach which is worth mentioning is based 
on a modified metric theory of gravity in which the usual 
Einstein-Hilbert action is replaced by a fourth-order ac- 
tion constructed with the conformal Weyl tensor [11, 12]. 
In particular, in the Newtonian limit, this theory leads to 
a gravitational potential in which, in addition to the 1 /r 
term, there arise two more terms proportional to r and 
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to r^. If the integration constants are suitably chosen, 
this potential leads to a satisfactory fit of the rotational 
curves of a large sample of galaxies without resorting to 
dark matter. 

However, the viability of the above alternative ap- 
proaches seems to be questioned by the relatively recent 
results of weak lensing observations of a unique cluster 
merger (the Bullet cluster), which point to the existence 
of dark matter independently of the assumptions regard- 
ing the nature of the gravitational force law [13]. 

Then, if we inspire ourselves with the Occam's razor 
principle, the failure so far to detect any of the hypo- 
thetical dark matter particles and the difficulties of the 
MOND-typc models lead us to pose the question whether 
there may be any chance that the missing mass prob- 
lem could be solved without introducing new particles 
or modifying general relativity. A possibility is that the 
missing mass is not actually composed of a non baryonic 
form of matter, but is rather a dynamical effect which 
occurs also in a vacuum, due to the local curvature of 
spacetime. Such a curvature would manifest itself as a 
gravitational force on test particles and light, mimicking 
the gravitational action of ordinary matter. 

There arc a few attempts pursuing this direction. For 
example, it has been shown in [14-16] that general rcla- 
tivistic nonlincaritics can play an important role even in 
the context of weak gravity, indicating, in particular in 
idealized models of galaxy clusters that the higher-than- 
expected constituent velocities can be accounted for by 
the baryonic mass alone without resorting to additional 
nonluminous matter. Similarly, it has been recently pro- 
posed that at least part of the discrepancy between the 
luminous matter and the total mass in a typical cluster 
can be explained by considering the relativistic influence 
of the nonhomogeneous distribution of far galaxies [17]. 
Another recent attempt [18] indicates the possibility that 
at least part of the action of dark matter could be ex- 
plained as a purely relativistic inertial effect due to the 
non-Euclidcan nature of three-space near ordinary mat- 
ter. Finally, it has been proposed in [19-23] that dark 
matter (and dark energy) may simply be a manifestation 
of the inhomogeneous nature of the cosmological energy- 
momentum tensor: in these papers, the authors show 
that, after an averaging process, a new source emerges in 
cosmology that mimics both dark components, depend- 
ing on the scale of interest. 

In the spirit of the above attempts, we study in this 
letter a toy model based on a particular vacuum so- 
lution of Einstein equations in which a hypothetical 
"Minkowskian" observer, experimenting with the motion 
of test particles in the asymptotically flat region, con- 
cludes that he is immersed in a continuous and homoge- 
neous (time-dependent) distribution of matter. Then, in 
spite of operating very sophisticated detectors in order to 
evidence the nature of these "unseen particles," he will 
be frustratingly unable to register any signal: he is im- 
mersed in a vacuum and the local "dark matter" effect 
that he feels is a purely dynamical one, solely due to the 
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FIG. 1: p-t section of the spacetime described by metric (1). 
The shaded, infinite cylinder represents the asymptotic re- 
gion, accessible to an hypothetical "Minkowskian" observer. 
In the bands above the singularity, dt is spacelike. 



curvature of spacetime. 

Since our model has several bizarre features, wc do 
not claim it provides a realistic answer to the origin of 
dark matter. It merely serves as an indication that dark 
matter effects may not necessarily be due to extra parti- 
cles: in the absence of corroborating signals one should 
be open to investigating alternative possibilities. 

We employ geometrized units in which c = G = 1. The 
Minkowski metric tensor is (1,-1,— 1,-1). 



II. MODEL 



Consider the following metric 
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where m>0, t>p>0, z & (—00,-1-00), G [0, 27r), 
and r ~ ^t'^ — ■ It can be obtained starting from the 
Schwarzschild metric 
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and applying the complex map 

Z ■■ (t, r, 6*, (j)] m) {it, ir, ix, 0; ~im) . (3) 
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This leads to the new metric 



A. Particle dynamics 



2m \ \ r 



1 



which, introducing the new coordinates 
t — r cosh X , p ~ r sinh x , z = t . 



(4) 



(5) 



takes the form (1). Metric (1), as evidenced when ex- 
pressed in the form (4) (with m G M), is the most gen- 
eral hyperbolically symmetric vacuum solution of Ein- 
stein equations. It was discovered by Harrison [24], and 
some of its properties have been investigated in [25]. Note 
that the sign of m in Eqs. (3)-(4) is the opposite of that 
considered in these references. Indeed, even though met- 
ric (4) is a solution of the vacuum Einstein equations for 
any value of m, the sign of the latter changes some fea- 
tures of the metric. For example, if m > there are no 
horizons, whereas if m < 0, a horizon at r = 2\m\ occurs. 
The nonvanishing curvature invariants of metric (1) are 
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showing it to be regular everywhere except at the 
singularity r = 0, precisely as it happens with the 
Schwarzschild one. However, whereas the Schwarzschild 
metric is static outside the horizon, metric (1) is nonsta- 
tionary in the in the region in which dt is timelike. 

Furthermore, whereas the Schwarzschild metric is 
asymptotically Minkowskian in space, metric (1) is 
asymptotically Minkowskian in time. Indeed, consider 
for the sake of definiteness an observer living in the re- 
gion where dt is timelike, and having access only to a 
radially limited portion of spacetime p < R, for example 
(see Fig. 1). Let (t, z, p, </)) be an event in such a region 
and, while keeping the coordinates z^p^(j) fixed, let the 
time t flow. Then, for large values of t/2m, (1) becomes 



dt^ 



dz^ 



{dp^ + p^c 



(8) 



In other words, whereas a local Schwarzschild observer 
becomes Minkowskian if he travels sufficiently far away 
in space (r — > oo), an observer located at some point 
in space in the manifold (1) becomes Minkowskian if he 
waits long enough (i — >■ oo, hence r — >■ oo since p has 
been fixed). 

Equation. (8) shows that the markers {t,z,p,(j)) can 
be interpreted as cylindrical coordinates for the asymp- 
totic or "Minkowskian" observer. By this interpretation, 
metric (1) displays an evident cylindrical symmetry: at 
any fixed time t, the space can be imagined as an infi- 
nite cylinder, having z as its axis, and which is radially 
limited by the singularity at Pmax = t. See Fig. 1. 



Consider an observer, who lives in the Minkowskian 
region. He probes the presence of matter in his neighbor- 
hood by observing the motion of test masses in free fall. 
Their trajectories are the timelike geodesies of metric (1). 
dz and are Killing vectors, so that the quantities 
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are conserved (the dot denotes the derivative with respect 
to proper time). Iz is the z component of the angular 
momentum per unit mass. The second equation shows 
that the velocity component z = u^t is conserved only 
in the asymptotically flat region. However, planar orbits 
with i = are always allowed, and these are precisely the 
type of motions we are interested in. Thus, a test particle 
in these conditions will orbit on a plane z — const, feeling 
an inward coordinate acceleration given by 



d^ 
dt^ 



+ K + 2Tly +r'v^v'^]vP , (10) 



In a "Newtonian" approximation, that is, for sufficiently 
low coordinate velocities [26] one has 



dp 
dt^ 



tt 



(11) 



where the Christoffel symbol T^f is expanded at the lead- 
ing order in 2m/ 1. Then, a test mass in planar motion 
feels the same force as if it were immersed in a cylindrical 
and homogeneous distribution of matter, which extends 
throughout the whole space, with a time-dependent den- 
sity 
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1 



G 27r(ci)3 



(12) 



where now we have explicitly introduced c and G. The 
mass parameter is the quantity mc^/G, whereas m is 
the geometrical mass. The observer concludes that some 
form of matter is homogeneously spread in the portion 
of space he is living in, even though no matter is present 
at all. Moreover, this matter density is decreasing with 
the third power of time, hence he also concludes that this 
distribution is dynamic, as if it were linearly expanding in 
every direction. Similar conclusions hold from the study 
of the trajectories of light rays. 

Of course, one might argue that similar considerations 
could in principle be applied to the Schwarzschild metric 
too, since the latter is also a vacuum solution. However, 
the difference between the two cases is crucial: in the 
Schwarzschild case the putative mass is concentrated in 
the origin, far away from the Minkowskian observer; in 
the case of metric (1), such "mass" is spread everywhere 
and the Minkowskian observer is immersed in it. 
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One can gain more physical insight about the features 
of the planar orbits as follows. The angular acceleration 
for a test-particle, under the same approximations which 
led to Eq. (11), is 







(13) 



then, ~ const. To be consistent with our approxima- 
tions, this constant has to be very small. Thus, the test 
particle rotates very slowly, with an almost vanishing an- 
gular velocity. Concerning the radial motion, the general 
solution of Eq. (11) is 



C.vlY.(ij), 



(14) 



where t = t/m and Ji and Yi are the Besscl functions 
of the first and second kind, respectively. Ci and C2 are 
integration constants, which are determined by the initial 
position po and radial velocity vq. Taking the derivative 
of Eq. (14) with respect to t we have, for t — > 00, 



vo 



mvQ 



O t 



-3/2 



(15) 



Equation (15) represents the asymptotic radial velocity 
of a test particle in our approximations. Then, a hypo- 
thetical probe in such a regime experiences an attractive, 
though rapidly decreasing, force, which lowers its radial 
velocity. 



III. DISCUSSION 

There are no indications that metric (1) has any di- 
rect physical significance. However, it often happens 
that even the study of nonphysical solutions may pro- 
vide insight on hidden aspects of a physical theory. In the 
context of general relativity, metric (1) is a clear exam- 
ple. Indeed, it allows for local, dynamical effects, which 
do not have an analogue in classical Newtonian mechan- 
ics: a Minkowskian observer, under rather general con- 
ditions, concludes that test particles are subject to an 
inward force which, according to Newtonian mechanics, 
he is led to interpret as due to a cylindrical and homo- 
geneous (though time-dependent) distribution of matter. 
However, metric (1) is a vacuum solution and, therefore, 
these effects are purely ascribable to the local curvature 
of spacetime. As such, this "dark matter" is fictitious 
and no experiment ever will be able to detect it. 
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FIG. 2: Planar timelike orbit [in (p, (j>) coordinates]. The test- 
particle trajectory is bent by the spacetime curvature even 
though the spacetime is empty. 



Appendix: Exact motion 

In order to obtain the effective distribution of matter of 
Eq. (12), we have imposed rather restrictive constraints, 
namely (i) the motion occurs in the quasi-Minkowskian 
region and (ii) the test-particle travels with low velocity. 
Indeed, if these requirements are not fulfilled, the motion 
could be quite different from the one caused by a cylin- 
drical and homogeneous distribution of matter. However, 
dynamical dark-matter-like effects still occur in a vac- 
uum, though they are not described by a simple effective 
potential such as the one we have considered in Eq. (11). 

Thus, in this appendix, we briefly discuss the exact 
motion of a test particle in metric (1). For this purpose, 
it is convenient to employ coordinates (r, r, x, (j)). In fact, 
in this case, the geodesic equations can be reduced to 
quadratures. For the case of planar motion (that is, z 
constant), the trajectory of the orbit is determined by 
the equations 



t>0 = h 



± 



dr 



dx 



sinh xVQ(x) 



(A.l) 
(A.2) 



where R{r) 
Z^/sinh^ X- K 
is analogous 
mass for the 
for timelike g 
A numerical 
a convenient 
displayed in 



:= r{r + 2m){K + m'^r'^), 6(x) := k — 
is an additional constant of motion, which 
the the total angular momentum per unit 
standard Schwarzschild geodesies; mo = 1 
leodesics, while mg = for the lightlike ones, 
example of the resulting timelike orbits, for 
choice of k, Iz and initial values of x, r, 0, is 
Fig. 2. 



5 



[1] J. Feng, Annual Review of Astronomy and Astrophysics 
48 (2010), in press, 1003.0904. 

[2] A. Giuliani, in Dark Matter and Dark Energy: a chal- 
lenge for Cosmology, edited by S. Matarrese, M. Colpi, 
V. Gorini, and U. Moschella (Springer, New York, 2010), 
to be published. 

[3] J. H. Oort, Bull. Astron. Inst. Netherlands 6, 249 (1932). 

[4] F. Zwicky, Helv. Phys. Acta 6, 110 (1933). 

[5] A. Dar, Astrophys. J. 449, 550 (1995), astro-ph/9504082. 

[6] C. Alcock et al. (MACHO), Astrophys. J. 542, 281 
(2000), astro-ph/0001272. 

[7] A. Masiero, in Dark Matter and Dark Energy: a chal- 
lenge for Cosmology, edited by S. Matarrese, M. Colpi, 
V. Gorini, and U. Moschella (Springer, New York, 2010), 
to be published. 

[8] M. Milgrom, Astrophys. J. 270, 365 (1983). 

[9] J. D. Bekenstein, Phys. Rev. D 70, 083509 (2004). 
[10] R. Reyes et al.. Nature 464, 256 (2010), 1003.2185. 
[11] P. D. Mannheim, Astrophys. J. 419, 150 (1993), hep- 
ph/9212304. 

[12] P. D. Mannheim, Prog. Part. Nucl. Phys. 56, 340 (2006), 

astro-ph/0505266. 
[13] D. Clowe et al., Astrophys. J. 648, L109 (2006), astro- 

ph/0608407. 

[14] F. I. Cooperstock and S. Tieu, Modern Physics Letters 
A 21, 2133 (2006). 



[15] F. I. Cooperstock and S. Tieu, Int. J. Mod. Phys. A22, 

2293 (2007), astro-ph/0610370. 
[16] F. I. Cooperstock and S. Tieu, Mod. Phys. Lett. A23, 

1745 (2008), 0712.0019. 
[17] A. Carati, S. L. Cacciatori, and L. Galgani, Europhys. 

Lett. 83, 59002 (2008), 0903.1355. 
[18] L. Lusanna, J. Phys. Conf. Ser. 222, 012016 (2010), 

0912.2935. 

[19] T. Buchert, Gen. Rel. Grav. 32, 105 (2000), gr- 
qc/9906015. 

[20] J. Larena, J.-M. Alimi, T. Buchert, M. Kunz, and P.-S. 

Corasaniti, Phys. Rev. D 79, 083011 (2009). 
[21] X. Roy and T. Buchert (2009), arXiv:0909.4155[gr-qc]. 
[22] A. Wiegand and T. Buchert (2010), arXiv: 1002.3912 

[astro-ph.CO]. 

[23] E. Kolb, S. Matarrese, and A. Riotto, New Journal of 

Physics 8, 322 (2006). 
[24] B. K. Harrison, Phys. Rev. 116, 1285 (1959). 
[25] M. Gaudin, V. Gorini, A. Kamenshchik, and 

U. Moschella, Int. J. Mod. Phys. D15, 1387 (2006), 

gr-qc/0512122. 

[26] In order for this approximation to be consistent, one must 
consider velocities such that v-' <^ 2m/t, as it easily fol- 
lows from an analysis of the various terms in Eq. (10). 



